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Abstract: A possible connection between the abundances of baryonic and dark matter 
^ (DM) has been explored so far mostly in the context of the so-called asymmetric DM. 

I— I Recently, a very different mechanism, dubbed "WIMPy baryogenesis" , has been proposed 

I to relate the baryon asymmetry to DM annihilation. The DM candidate is a weakly 

interacting massive particle (WIMP), and the usual WIMP scenario is slightly extended to 
accommodate baryogenesis, which is accomplished around the time of DM freeze-out. We 

^> construct an effective field theory that encompasses a quite general class of models which 

implement the WIMPy baryogenesis. Under some reasonable, simplifying assumptions, 
we show that a good portion of the parameter space is allowed for these models, after 

CN experimental constraints are taken into account. Bounds from the LHC require that the 

WIMP be heavier than 400 GeV. 
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1 Introduction 

The presence of non-luminous and non-baryonic matter, the so-called Dark Matter (DM) [1- 
7], and the existence of the baryon asymmetry of the Universe (BAU) [8-10] are two well 
established facts. Cosmic microwave background (CMB) anisotropy observations by the 
Wilkinson Microwave Anisotropy Probe (WMAP) yield an accurate determination of the 
total amount of baryonic matter [9, 10], 

Ob/i^ = 0.02260 ± 0.00053 , (1.1) 

and of non-baryonic matter 

f^D^/i^ = 0.1123 ± 0.0035 . (1.2) 

The fact that the two abundances are comparable 0dm /^b ~ 5 can be fortuitous, or may 
be the sign that they have a common origin. This intriguing possibility has been vastly 
explored in the literature, mostly by invoking asymmetric DM (ADM) scenarios [11-47]. 
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The common feature of all such models is that the DM abundance is determined by a 
matter- antimatter asymmetry in the dark sector, which in turn is connected to the baryon 
asymmetry in the visible sector. The DM asymmetry can be produced prior to the BAU, 
can emerge from it, or can be produced during parallel and competitive processes, but in all 
cases the observed non-baryonic matter results from an asymmetry, implying the necessary 
suppression of the symmetric component. The ADM scenarios require a non-trivial dark 
sector. An argument in support of this choice is that the visible sector is rich and complex, 
why should the dark one be much simpler? Although this logic is perfectly sensible, we must 
admit that we still know very little about the nature of DM. Therefore it seems reasonable 
to first tackle the problem keeping the dark sector as minimal as possible. A very simple 
framework is that where the DM candidate is given by just one weakly interacting massive 
particle (WIMP). Motivated by the well-known "WIMP miracle", one can then ask the 
following question: is it possible to extend minimally the WIMP scenario in order to make 
a connection between the DM and baryon abundances? 

A few attempts in relating the relic density of such symmetric WIMP to the observed 
BAU have recently been made. The authors of [48] proposed a mechanism that they 
dubbed "dark matter assimilation". The idea is that singlet DM, (they study a bino-like 
neutralino as an example) is efficiently depleted in the early Universe by being assimilated 
into new quasi-stable heavy states, ^ and ^ via the reaction ~^ ^'P^ where (/> is a 
Standard Model (SM) particle. These new heavy states carry the baryon asymmetry. 
Their subsequent annihilations and decays (into DM and SM quarks) yield the correct 
DM relic density and BAU. The minimal version of such models is quite economical as it 
requires only the addition of singlet DM and two new heavy states on top of the SM particle 
content. Another mechanism, that goes under the name of "baryomorphosis" and shares 
some similarities with the DM assimilation, was studied in [26, 38]. In the baryomorphosis 
scenario a large asymmetry is originally stored in a heavy scalar field, which then decays 
into two colored scalars, the annihilons, with mass at the TeV scale. Their late annihilations 
into singlet scalar DM particles set the correct DM abundance and baryon asymmetry. 

More recently, Cui, Randall and Shuve introduced a new mechanism, which they refer 
to as "WIMPy baryogenesis" [49]. Here the DM is a Dirac fermion^ with a mass of at least a 
few hundreds GeV. One also needs to add two new states, one heavy (~ TeV) charged under 
the SM quantum numbers, ip, the other one very light and uncharged, n. DM annihilates 
into if) and a quark (or a lepton), and ifj subsequently decays into the light, sterile state n, 
storing the negative asymmetry in baryon (lepton) number in a sequestered sector. The 
DM relic density is that of a thermal WIMP, and by the time the DM annihilations freeze 
out one has generated a baryon asymmetry, either directly, when the annihilation is into 
■0 plus quark, or via leptogenesis, when the annihilation is into if^ plus lepton. In the latter 
case the lepton asymmetry has to be generated before the electroweak phase transition 
(EWPT) so that sphaleron processes can transfer the asymmetry to the baryon sector, 
with the consequence that the masses of DM and -0 must be at least O (TeV) . In the case 

^This is the possibility studied in their explicit models. A Majorana fermion or a scalar can in principle 
work as well. 
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where DM annihilates into quarks, baryogenesis can occur after the EWPT and masses 
down to a few hundreds GeV are ahowed. 

There is a conceptual difference between the WIMPy baryogenesis framework and the 
typical ADM frameworks. Whereas both aim at predicting the correct baryon and DM 
abundances, they stem from different motivations. ADM models are built in order to ex- 
plain the observed ratio Q^M/^h ~ 5, understood as non- fortuitous. WIMPy baryogenesis 
models, instead, are based on taking seriously another coincidence: the WIMP miracle. 
Adding only a few new ingredients, one keeps such a miracle and can easily accommodate 
baryogenesis. On the one hand, ADM gives up the WIMP miracle; on the other hand, 
WIMPy baryogenesis does not explain why Qdm/^Y) ~ 5, but as we show in this work, for 
values of the parameters that we consider natural, and are still well within experimental 
constraints, these models give the right numbers for the DM abundance and the BAU. 

The aim of the current paper is to study a general class of models that implements the 
WIMPy baryogenesis mechanism. In the spirit of keeping the models as simple as possible 
we identify the minimal particle content that does the job. In Ref. [49] , the possibility of 
DM annihilations generating a lepton asymmetry, which is then converted into a baryon 
asymmetry through sphalerons, was somewhat emphasized. In this work we only consider 
models where DM annihilates into a quark plus an exotic, heavy antiquark, thus producing 
the BAU directly. We find the latter scenario more appealing than the former for a few 
reasons. First, for the sake of simplicity, we can avoid the extra step with sphalerons needed 
in the leptogenesis case. Second, because the exotic quark is colored, as the name suggests, 
the LHC has a chance of discovering it or, alternatively, it could put some severe bounds 
on the model. Third, given that the DM directly couples to quarks, we can in principle 
hope for a direct detection signal, which would certainly not be there if DM interacted 
only with leptons. Thus, the scenario we study here seems more testable than the WIMPy 
leptogenesis, although the latter still remains a possibility worth exploring. 

Our approach is a little different than the one in [49]. They build a UV-complete 
model, including pseudoscalars that mediate the DM annihilation and can have masses 
of the same order as the DM and the exotic quark; the Lagrangian they use as the basis 
for their calculations only includes renormalizable interactions. We do not include any 
mediators in our model instead, assuming that they are all much heavier than the fermions, 
and in the spirit of an effective field theory (EFT) we write only four-fermion interaction 
terms. Doing so we have a more generic and richer class of models that allows for new 
operators and for new DM annihilation channels. 

The paper is organized as follows. In Section 2 we present the particle content and 
the Lagrangian of the model, then we discuss how the baryon asymmetry is generated. In 
Sections 3, 4 and 5 we study the experimental constraints on the model from the LHC, 
cosmology, and DM direct detection. We summarize our conclusion in Section 6. Some 
technical details are included in the final appendices. 
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2 An effective-field-theory approach 

2.1 Field content of the model and Lagrangian 

We add to the SM the minimal particle content which is needed in order to have a successful 
WIMPy baryogenesis. All the new particles are fermions. We consider vector-like gauge 
singlet DM X and X, vector-like exotic quark color triplets ip and ■i/', and a massless singlet 
n, into which the exotic quark decays, as we will explain in the next subsection. 
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Table 1. Particle content of the model, u and d are the right-handed up and down quarks of 
the SM. The rest of the SM quarks also have charge —1, while all the leptons are neutral under the 
Z4 symmetry. The reason for these charge assignments is explained in Appendix B. 

Notation and conventions We use the two-component spinor formalism and we follow 
closely the conventions of Ref. [50]. The advantage of such a formalism versus the four- 
component spinor one is that Fierz identities are easier, which greatly simplifies the task 
of finding a complete, linearly independent basis of dimension six operators. Occasionally, 
it proves convenient to switch back to four-component notation, in which case we denote 
the spinors as follows 

From eq. (2.1) it is evident that X and X^ represent DM particles, while X and X^ anti- 
DM particles. 

A discrete symmetry is needed in our model, in order to stabilize the DM and to avoid 
dangerous decays of the exotic quark that could spoil the baryon asymmetry. The reason for 
the choice of the discrete charges assigned to the various fields is discussed in Appendix B. 
In the same Appendix we also show why, restricting ourselves to global discrete Abelian 
groups, the minimal choice within our field content is Z4. We emphasize the fact that the 
Z4 we impose is not generic for the WIMPy baryogenesis mechanism. Rather, it is tied 
to our specific models. It is possible that other models with different symmetry groups do 
the job, although we have not investigated such a possibility. 

In order to study the phenomenological implications of this model, we write down 
an effective Lagrangian that includes all the dimension six operators Oi (four-fermion 



-4- 



operators) consistent with the field content and the quantum numbers hsted in Table 1, 

£d^5;A?0.. (2.2) 

i 

Here we have chosen to parametrize the couplings as (Aj/A)^, with dimensionless Aj's and 
a fixed mass scale A. The list of 20 (plus Hermitian conjugates) operators Oi is given in 
Appendix A. Given the particle content we consider, this set of operators, seemingly just 
from scalar exchanges, spans a complete and irreducible basis. All the possible vector- 
(e.g. {X'^ a^^ X){ua fj^'ip'^)) and tensor- (e.g. {Xa^^^ X){'il)a^j_uu)) operators can be related to 
the scalar ones in our list by Fierz identities. 

The EFT approach is valid only as long as the biggest momentum involved in the 
processes we are considering, fcmax, is such that /cmaxx < 1- In our study, A^max is given by 
some temperature, T, before the DM freezes out. As we will explain in the next subsection, 
the baryon asymmetry should build up between the time the washout processes freeze out, 
z = zyfo, and the time of DM annihilation freeze-out, z = zpo, where z = m^/T. As in 
usual thermal WIMP scenarios, z-po ~ 25, while we find that typically 2;wo ~ 10. Given 
that we scan over a DM mass range up to 1 TeV, with zwo = 10 giving the reference value 
for the highest temperature, we have /cmax ~ ""t-xZ-^^wo ~ 100 GeV. The condition for the 
validity of the EFT approach then translates into a bound for the couplings, namely 

^ < (100 GeV)-i . (2.3) 

In the following numerical evaluations we fix A = 10 TeV, which in turn implies Aj < 100 ^. 

To keep the number of parameters in the numerical analysis manageable, we set some 
equalities among the relevant couplings and we relabel them for the ease of the discussion'^: 

• \si = Ai = A3, As2 = A2 = A4: couplings for DM (or anti-DM) annihilation into 
ipu and ifj'^v) in the s-channel; 

• At = A5 = Aq: couplings for DM (or anti-DM) annihilation into and ijj'^v) in the 
t-channel; 

• Awo = Ag = An: couplings responsible for the "pure" washout processes^; 

• Xn = A20: coupling implying the decay ijj — )• ddn. 

^ This last bound can lead to some confusion, thus a comment is in due order. In the parametrization 
we chose for the Lagrangian (2.2), the couplings Ai can be thought of as dimensionless Yukawa's. In a given 
UV-completion of our theory, they would be subject to the usual perturbative bound, Ai < An, which seems 
to be in conflict with our Ai < 100. The point is that the only sensible condition for our EFT is expressed 
as ^ < (100 GeV)~^. One could always keep Ai below 47r by lowering the scale A. In other words, the 
bound Ai < 100 is just a consequence of our parametrization and the choice of fixing A to 10 TeV. Values 
of Ai bigger than An can just be thought of as lowering the scale A below 10 TeV. Different parametrization 
would obviously lead to the same physical results. For example, other DM studies based on EFT approach 
(e.g. [51-53]) assume instead 0(1) coefficients and constrain one scale Ai per operator. 

^In Appendix A the reader can find the operator Oi associated to Ai, with i = 1 . . . 20. 

*The "pure" and "mixed" washout processes are defined in the next subsection. 
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The couplings A^i, Ac,2 and At are also involved in the "mixed" washout processes, that are 
related to the annihilations by crossing symmetry. The reason for keeping two couplings in 
the s-channel is evident if we write the corresponding operators in four-component notation 

^ [(A^i + X's2){xx'' + rx) + (A^i - A22)(x75X^ - ri5X)] [UPl^] • (2-4) 

This expression shows that, for general A^i and As2, we have a mixture of scalar and 
pseudoscalar channels. In our analysis we will consider two limiting cases: 

As = A^i = +Ag2 scalar s-channel, (2-5) 
Ap = A^i = — Ag2 pseudoscalar s-channel. (2-6) 

It is also worth noticing here that what we define as t-channel is equivalent to a vector 
exchange in the s-channel. This can be understood simply by using a Fierz identity: 

A? 



lA? 



= "2^ + {xI^Plx'WPlI^'^)] . (2.7) 

In the last line, we have translated into four-component notation. 

Many couplings do not appear in the relabeling list above. Two of them, A7 and Ag, 
contribute to DM annihilation into two quarks and to DM direct detection, and will be 
discussed in Section 5. The couplings for the annihilation XX — )• '\\)'^ 

A^^ = Ai3 = Ai4 = Ai5 = A16 = Ai7 = A18 = Ai9 (2.8) 

are assumed to be small and are neglected in this study. This channel, corresponding 
to a CP-conserving DM annihilation and kinematically open when m-^ > m^, does not 
contribute to the BAU generation and reduces the efficiency of the WIMPy baryogenesis 
mechanism, so it has to be suppressed compared to the other annihilation channels. We 
require A^^ ^ Asi,2) Aj. The remaining two couplings, Aio and A12, are included in the list 
in Appendix A for completeness, but do not play any role in the rest of our discussion. 

2.2 Generation of the baryon asymmetry via WIMP annihilations 

DM annihilates into a u quark and an exotic quark, ip. We emphasize that in these 
models the DM, despite being a Dirac fermion, annihilates with itself rather than with 
its antiparticle. This can be seen explicitly from eq. (2.4), for example. During these 
annihilations, CP is violated and an u-quark number asymmetry is created, in an equal 
and opposite amount to a ^-number asymmetry^. After being produced, the exotic quark 



^In principle, one could worry about sphaleron processes even in this context. If the BAU generation 
occurs before EWPT, the baryon asymmetry is converted into a conserved B-L charge, while the B charge 
is set to zero by sphalerons. After their freeze-out, the B-L number is split into B and L number. However, 
given the masses considered in this work, the baryon asymmetry generation occurs after EWPT, with no 
reprocessing effects from sphalerons. 
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decays into a SM singlet fermion, n, plus two antiquarks: ip — >• ddn. It is crucial that ip 
does not decay only into SM quarks, since that could eliminate the asymmetry. The Z4 
symmetry prevents such dangerous decays. To avoid overdosing the Universe, n has to be 
light^. For the sake of simplicity we take it to be massless. Two scenarios are contemplated: 

1. The singlet n carries baryon number +1. In this case the decay ip — >• ddn is baryon- 
number-conserving, but n is sequestered in a sterile sector, so we are left with a net 
baryon number in the visible sector. The overall process violates the SM baryon 
number. 

2. The singlet n does not carry baryon number. The decay of "0 explicitly violates 
baryon number and it contributes to the baryon asymmetry. 

Both cases satisfy the first Sakharov condition [54], that is baryon number violation. The 
other two conditions are also easily satisfied: CP violation is achieved with complex cou- 
plings Aj and with the interference between tree-level and one-loop diagrams (see Ap- 
pendix D); departure from thermal equilibrium is automatically implemented, given that 
WIMP annihilation around freeze-out is out of equilibrium. 

There are processes that can potentially wash out the asymmetry. They are shown 
schematically in figure 1. The processes ijju — t- XX and uX — t- X'^ip'^, involving DM 




Figure 1. Washout processes. We show mixed washout processes in the upper two diagrams 
and pure washout processes in the lower two. 

particles, are referred to as "mixed" washout. They are obtained from the DM annihilation 
diagrams by crossing symmetry, and so involve the same couplings A^i, As2, Af. The 
processes in the two lower diagrams, that involve only ip and u, go under the name of 
"pure" washout. We will discuss further about these processes in Section 4. 

®An interesting possibility is that n constitutes the extra degree of radiation at BBN. We will explore 
this in more detail in future work. 
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A very important result, emphasized by the authors of [49], is that in order to produce 
a significant baryon asymmetry, washout processes must freeze out before WIMP freeze- 
out. To achieve this early washout freeze-out, one needs either a ip heavier than the DM, 

> m-^, so that the washout is Boltzmann suppressed while DM is still annihilating, or a 
small couplings, such that the washout cross section is small compared to the annihilation 
cross section. 

3 Constraints from the LHC 

One of the new particles that we need in our models, ij:, is colored, which makes it a good 
candidate to be discovered at the LHC. Alternatively, the LHC can put severe bounds on 
these models. ip can be pair-produced at the hadron collider, pp — )• ipip. Given that each 
-0 decays into two (anti-)down quarks plus a singlet n, the signature to look for is four jets 
plus missing energy. Both the CMS [55] and the ATLAS [56] collaborations search for such 
a signature in the context of supersymmetry (SUSY), and they put bounds on the masses 
of gluinos and squarks from the process pp ^ gg ^ 4j -1- where the missing energy 
is carried away by the lightest neutralinos. In the case of a massless neutralino, which is 
considered in [56], the SUSY process is completely analogous to the one we are interested 
in, with tp in place of the gluino. The production cross section differs only by a group 
theory factor, since ■i/' is a color triplet, while the gluino is an octet. Once that is taken 
into account, we can translate the bound^ on the gluinos, rUg > 960 GeV, into 

> 800 GeV. (3.1) 

In these models DM has to annihilate into a quark, whose mass can be neglected, and 
an exotic antiquark, tp, in order to generate a baryon asymmetry. Therefore we have the 
kinematical bound 

< 2 m^. (3.2) 
Thus, the DM has to be heavier than ~ 400 GeV. 

4 Constraints from cosmology 

In this section we assess the impact of the DM relic density and BAU constraints on 
our parameter space. To that end, we make use of Boltzmann equations, detailed in 
Appendix C, for the evolution of the DM density and the baryon asymmetry. 

In the numerical evaluations, we fix A = 10 TeV, and we scan over and m^, with 
< 1 TeV, and 800 GeV < < 2 m^^. The lower and upper bounds on are 
dictated by the LHC and the kinematics of DM annihilations, respectively. We recall that 
the dimensionless couplings Aj's have an upper bound of ~ 100, as required by the validity 
of the EFT approach. 

^This is the bound from [56] when the squark is very heavy, which is appropriate in our EFT context, 
where all the mediators are very heavy and have been integrated out. 
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4.1 Dark matter relic abundance 



The Boltzmann equation governing DM annihilations is given in eq. (C.l). As we work 
in a first order expansion in the asymmetries, the small contributions proportional to e.g. 
e X Yau are neglected, and the DM relic density is set by the typical thermal freeze-out. 
In principle the annihilations can have a CP-conserving (CPC) and a CP-violating (CPV) 
contribution. The CPC channels, governed by the couplings A^^, A7 and Ag reduce the 
WlMPy baryogenesis efficiency, so they have to be subdominant if one wants to generate 
a sizable baryon asymmetry. On the one hand, A7 and Ag can be somewhat suppressed 
to avoid conflict with direct detection searches, as we explain in section 5. On the other 
hand, we have to assume A^^ <C Asi^2, Aj when this channel is kinematically open, i.e. for 
m^/m^ < 1. 

It is instructive to inspect the CPV DM annihilation rate^. In the low temperature 
limit (large z = m^/T) it reads 



where x = m^/(2m^). Eq. (4.1) shows the known result that the scalar s-channel is 
velocity suppressed compared to the pseudoscalar one. Since in the rest of this section 
we will analyze the scalar, pseudoscalar, and t-channel cases separately, it is worth having 
in mind their relative contributions. To this end, we display in figure 2 the rate 
(normalized to n^{z)H{z)) for various limiting cases. It is clear that the pseudoscalar and 
t channels dominate. Thus, a larger As, compared to Ap and At, will be needed to get the 
right relic density. The kinematic suppression of the rates is also manifest in figure 2 for 
m^/m^ ^1-8: the larger the ratio m^/m^ is, the larger the couplings will have to be in 
order to compensate for such a suppression. 

The expected constraints are confirmed in figure 3, which shows contour levels for 
the coupling Ap (upper left pane), A^ (upper right pane) and A^ (lower pane) needed for 
generating the DM relic density abundance measured by WMAP, in the [m^/rriy., m^] 
plane. The white lower left region in each plot corresponds to a i/' mass already ruled out 
by the LHC (m^ < 800 GeV). In general, the couplings Ap^^^t have to be larger than 10. 
We see that for a fixed m^/m^ ratio, smaller couplings are needed for larger m-)^. Indeed, 
with the increase of the DM relic abundance gets reduced due to the thermal average 
(for a fixed DM annihilation cross section), and to the increase of the DM cross section 
itself, hence smaller couplings are required. Oppositely, for a fixed DM mass, as shown 
in figure 2, the annihilation cross section decreases for larger ^0 masses, therefore larger 
couplings are needed in order to compensate for the reduction of the available phase space. 

Combining pseudoscalar and t-channel (Ap = Xt and Ag = 0) or scalar and t-channel 
(As = Xt and Ap = 0), we find the pattern that would be expected and inferred from the 
limiting cases analyzed above: the pseudoscalar dominates over the t-channel, and the 
t-channel dominates over the scalar one. This is illustrated in figure 4. 

*See Appendix C for the definition of the rate. 





3 



(4.1) 
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m, =800 GeV, y\=10 TeV , z=15 

50r — ■ — ' — ■ — • — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ' — ■ — — 




L h ........... . . , . . . v1 

1.0 1.2 1.4 1.6 1.8 2.0 



Figure 2. Dark Matter annihilation rate as a function of the mass ratio m^/m~^, for 
various limit cases. From top to bottom, the curves represent the pseudoscalar s-channel, the 
t-channel and the scalar s-channel. For illustration, all respective couplings have been set to 1. 

The results in this subsection are in Hue with those obtained in other EFT ap- 
proaches [51-53]. 

4.2 Baryon asymmetry 

Now we want to impose the constraint from the BAU, on top of the DM rehc abundance 
already considered in the previous subsection. The final baryon asymmetry, that has to 
match the measured value, results from competing processes: CP-violating DM annihila- 
tions generate an asymmetry, while washout processes tend to deplete it. 

Let us first consider the asymmetry, e, defined in eq. (D.l) as the difference between 
the rate for DM annihilation into t/^u and the rate for the annihilation into ip^u^ , normalized 
to the sum of the rates. As we explain in Appendix D, after we define the couplings As, Xp 
and Xt by setting equalities among the a priori different Ai's, we need to make the choice 
of assigning the CP-violating phases to some of the couplings. The most economical choice 
is to assign just one phase to the pure washout coupling, Awo^ 

Awo = |Awo|e'^ = A9 = Aii. (4.2) 

With this choice, we have (see eq. (D.4)) 

eoc |AwoPsin(2 5). (4.3) 

To further simplify the analysis, we make the assumption that the phase is large and 
we set 6 = tt/A. This is not in conflict with current electric dipole moments (EDMs) 
measurements that would be the most constraining for these phases, given that the lowest 
order contribution to EDMs in these models is at three loops [49] and so very suppressed^. 

® Although the diagrams in our EFT contributing to EDMs are difTerent than the ones in [49], they are 
stiU three-loop suppressed. 
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Figure 3. Contour level for the couplings Xp (upper left pane), At (upper right pane) 
and As (lower pane) needed for generating the DM relic density in the [m^/m-^, m-^] 
plane. The white lower left part corresponds to < 800 GeV, which is excluded by the LHC. 



In figure 5 we show e for various limiting cases. Note that the pseudoscalar and scalar 
channels give higher values of e than the t-channel. In the cases with Xt = 0, the dependence 
on Xs and Xp cancels out in eq. (D.4), so that the scalar and the pseudoscalar s-channel 
asymmetries are equal. 

Next, let us consider the washout processes. The same couplings A^, Xp and At, that 
enter in the annihilations and have already been constrained by requiring the correct DM 
relic density, also contribute to the mixed washout. In figure 6 we show how the rates 
for the mixed washout compare in the three limiting cases (again normalized to 



n^{z)H{z)). Besides the very strong suppression of these rates with increasing m^/ 
from kinematic closure, it is interesting to note that the pseudoscalar case (blue line) gives 
the least washout. Combined with the high annihilation rate and with the high values 
e, as shown in figures 2 and 5, this makes such channel the most promising for achieving 
WIMPy baryogenesis. The pure washout rate ")^q is also depicted in figure 6, where we 
set all the couplings equal to 1 for illustration and comparison. It is important to keep in 
mind that each washout rate scales as the corresponding coupling Aj to the fourth power. 
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Figure 4. Contour level for the couplings \p — At (left) and Xg — At (right) needed for 
generating the DM relic density in the [mjj^/m-^, m~^] plane. The white lower left part 
corresponds to m^p < 800 GeV, which is excluded by the LHC. 



m^=800 GeV, A=10 TeV , z=15 

10-3 ^ 




1.0 1.2 1.4 1.6 1.8 2.0 



Figure 5. Comparing e's for various channels. The green dashed line is for As = 1 and 

Ap = A( = 0, which equals the case where Ap = 1 and A^ = Af = depicted with the blue solid line. 
The dotted red line is for At = 1 and A^ = Ap = 0, corresponding to a pure t-channel, but also to 
a t-channel plus scalar s-channel. The black long-dashed line stands for the case Xp — Xt ^ 1 and 
A, =0. 

We have already seen that the DM rehc abundance requires the coupUngs As, Ap and Xt to 
be at least of order 10. In figure 7 we see that AwOi instead, is typically of order 1. As a 
consequence, the mixed washout rates are enhanced by a factor of 10^ compared to pure 
washout rates in our processes, thus they are dominant. 

In figure 7 we depict contour levels for the modulus of the coupling Awo needed 
for generating the measured value for the BAU, in the plane [m^/m^, rriy.], again for the 
limiting cases considered before. Although not shown in the plots, the values for Ap, A^ and 



- 12 - 



m^=800 GeV, A=10 TeV , z=15 




Figure 6. Washout rates as a function of 111,^/171-^. The rates are normalized to n'^ (z) H (z) . 
The color code is the same as in figure 2: blue is for the pseudoscalar s-channel, red for the t- 
channel and green for the scalar s-channel. The black line here stands for the pure washout rate. 
For illustration, all respective couplings have been set to 1. 

At are fixed to reproduce the correct DM relic density. It is worth noticing that, except for 
the white region already ruled out by the LHC, the rest of the parameter space could give 
rise to the correct DM relic density and the baryon asymmetry, for reasonable values for the 
couplings. The behavior of |Awo| can be understood as follows. Recall that e oc |AwoP"^x' 
while the mixed washout, which dominates over the pure one, is proportional to A^ g For 
a fixed value of m^, increasing m^, both e and the washout decrease. But while the washout 
rates decrease quickly, e goes down slowly for m^/m^ < 1.8, so |Awo| has to decrease in 
this direction in order to not overproduce the asymmetry. For m^/m^ ^1-8 the washout 
processes are not important anymore and e would become too small, thus |Awo| has to 
invert the trend and start increasing. For a fixed m^/my. ratio smaller values of |Awo| are 
needed when increases. 

Figure 8 shows contour levels for e (in the low temperature limit) generated when 
imposing both the DM relic density and BAU constraints. The parameter e roughly follows 
the same behavior as if we were keeping Awo constant: it decreases with the increase of 
m^/m^ (see equation (D.4)). 

For completeness, we display the remaining two limiting cases, where Xp = Aj with 
As = 0, and As = At with Ap = in figure 9. The constraints on |Awo| and on e are 
shown in the upper and lower panels respectively. Whereas we had a clear dominance of 
the pseudoscalar channel in the determination of the DM abundance, we see in the figures 
of this subsection that the parameters e and Awo do not vary significantly among the three 
limiting cases. 
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Figure 7. Contour levels for the modulus of the coupling \wo needed for generating 
the measured BAU, in the [m^/m-^, m^^] plane. We display the pseudoscalar (upper left 
pane), t-channel (upper right pane) and scalar (lower pane) cases. 



5 Constraints from dark matter direct detection 



We have two operators in our effective Lagrangian that contribute to the direct detection 
of DM at tree level: 



1 

A2 



+ h.c. . 



When eq. (5.1) is translated into four-component-spinor language 



Ag - A 



4A2 



4A2 



[{X1^75X){UJ^,U) + (x7^75X)(t/7^75f/)] 



(5.1) 



(5.2) 



we see clearly that the first term (vector coupling) gives a spin-independent (SI) contribu- 
tion, the last term (axial-vector coupling) a spin-dependent (SD) contribution, while the 
second and third terms are velocity suppressed and can be neglected. Note that the SI 
term is proportional to the difference of the couplings, while the SD one to the sum. 
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Values of e for Xs=h=0 Values of 6 for As=Ap=0 




0.8L . . . . . J 

400 500 600 700 800 900 1000 
m^. [GeV] 

Figure 8. Contour levels for e. We display the pseudoscalar (upper left pane), t-channel (upper 
right pane) and scalar (lower pane) cases. 



The operators in eq. (5.1) are also responsible for DM annihilation into a quark plus an 
anti-quark. This annihilation channel, which does not contribute to the asymmetry, would 
be competing with the one into quark plus exotic anti-quark. We want the former to be 
suppressed with respect to the latter, in order to generate the correct BAU. Therefore, 
even strict bounds on the couplings A7 and Ag from direct detection, would not challenge 
these models. Put another way, for the WIMPy baryogenesis to work, A7 and Ag must be 
suppressed with respect to A^, Ap and At, which would explain why no signal has been seen 
so far in direct detection experiments. 

Let us take a look at the bounds, starting with the case that A7 and Ag are different 
and consider the SI contribution. The WIMP-nucleon cross section is given by 

do = , [0.6) 

TT 



In equation (5.3), B^^ = {A + Z) + aa {2 A — Z), A is the atomic mass number, Z is 
the atomic number and u^v is the WIMP-nucleus reduced mass. and ad correspond to 
the couplings between DM and the up- and down-quark respectively: a„ = -^j^, arf = 0. 
Figure 10 depicts the impact of the XENONIOO exclusion limits, after 225 live days of 
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Figure 9. Contour levels for the modulus of the coupling Xwo (top panels) and the 
parameter e (bottom panels) needed for generating both the DM relic abundance and 

the BA U. Left panels show Xp = At with As = 0, while on the right we set Xg = A* with Xp = 0. 



data [57], on the plane [A| — Ay, rriy.]. The upper orange region is ruled out. For the 
relevant mass region of the present study (i.e. > 400 GeV) the combination of the 
couplings Xj — Xg has to be smaller than ~ 10. 

If Ay = Ag, there is no SI contribution, but the SD one is at its maximum. The SD 
limits are a few orders of magnitude weaker than SI limits, and are of no interest for the 
present study, given that they are weaker than the bound from the validity of the EFT 
approach, Aj < 100. 

Constraints on the couplings A^, Ap and A^ would be more interesting, because they 
play a crucial role in the generation of the baryon asymmetry and the determination of the 
DM relic abundance. The lowest order contribution to the direct detection involving these 
couplings is naively expected to be at one loop. 

Let us take a closer look at this statement and let us consider first the situation where 
only the s-channel is turned on (A* = 0). There are two one-loop diagrams, as shown 

^°This is certainly true for lieavy DIvI, wliicli is our case, where the best hmits are from direct detection 
experiments, see e.g. [58, 59]. For light DIVI one would get better SD limits from colliders [60-62]. 
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Figure 10. XENONIOO sensitivity. The upper orange region on the parameter space [A| — 
Ay, m^] is already ruled out by the latest XENONIOO [57] measurements. 

in Figure 11. We show the As couphng in the figure, but if we replaced it with Ap the 
following argument would still apply. It is easy to see that, in the limit of zero external 




Figure 11. Diagrams for direct detection at one loop. In the limit of zero external momenta 
the sum of these two diagrams vanishes. 

momenta, appropriate for direct detection, the sum of the two diagrams vanishes. This is 
understood just by looking at the fermionic propagator for ip in the loop, the direction of 
which in the second diagram is opposite to the one in the first, thus giving a relative minus 
sign. This is enough to conclude that the contribution to the cross section is not only 
loop suppressed but also velocity suppressed. As a consequence current direct detection 
experiments place virtually no limits on the coupling A^. The exact same conclusion holds 
for the pseudoscalar coupling Ap. 

For the t-channel {Xt 7^ 0, = Ap = 0), there is no analogous obvious cancellation 
at the level of one- loop Feynman diagrams, but it turns out that there is only a SD 
contribution, which is loop suppressed. Since, as we said, the SD limits are not even 
important for the tree-level couplings, we are definitely far from putting constraints on At 
with direct detection. 
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6 Summary and discussion 



In this work we have investigated whether a general class of WIMPy baryogenesis models 
is viable, after experimental constraints are taken into account. Our models are based 
on the same mechanism and on the same external particles as in Ref. [49]. However, by 
following an EFT approach and writing down a complete list of four-fermion operators, we 
extend and generalize their study, considering all the possible DM annihilation channels. 
The models considered here require the presence of a heavy fermion, -0, which is crucial to 
the success of the whole mechanism. Because ip is colored, the LHC represents an excellent 
laboratory for testing these models. Although we have not yet studied in detail possible 
collider signals, current LHC searches already put a lower bound of 800 GeV on the mass 
of tp, which in turn directly translates into a lower bound of 400 GeV on the DM mass. 
With the impressive pace at which the LHC and the ATLAS and CMS collaborations are 
operating, this bound can increase relative quickly, pointing to even higher masses, or, in 
a better (luckier) scenario, a heavy colored fermion could be discovered soon, which would 
provide a hint that these models could be realized in nature indeed. 

In this work we focused mainly on the cosmological aspects and we examined in some 
detail the constraints from the measured DM relic density and BAU. We considered three 
different channels for DM annihilation into a quark and an exotic antiquark: scalar and 
pseudoscalar s-channel, and t-channel. We found the pseudoscalar channel to be the most 
promising: it has the highest annihilation cross section, the lowest washout cross section 
and it generates a large asymmetry e. This combination results in lower values of the 
coupling \p, compared to As and \t, and in the most efficient production of the BAU. In 
the spirit that lower rather than higher values of the couplings are generally preferred in 
the EFT, our analysis, in all cases, points toward a high DM mass, between 800 GeV and 
1 TeV, and a small hierarchy between ^ and Xi < 1.4 (see the figures in Section 4). 

We also considered bounds from direct detection. These constrain only two operators 
that would be responsible for the annihilation of DM into a pair of quarks. Given that this 
channel does not contribute to the generation of the baryon asymmetry, we want it to be 
suppressed anyway. In this sense such bounds do not challenge these models at all. There 
are in principle one-loop diagrams, involving the couplings that also enter the generation 
of the asymmetry, that could contribute to the direct detection cross section. We showed 
that they are not only loop suppressed, but also velocity suppressed. Thus, this scenario 
is out of reach for current direct detection experiments. 

Since CP violation is a crucial ingredient in these models, one has to worry that the 
physical phases are not too constrained. As already pointed out in [49], it seems that the 
strongest constraints on the phases would come from EDMs measurements. For the models 
of [49], the lowest order contribution to the neutron EDM only appears at three loops, and 
as a consequence their phases are not much constrained at all. In our EFT context, the 
diagrams contributing to the neutron EDM are slightly different, but it is still true that the 
lowest order contribution appears at three loops, so we reach the same conclusion: the CP 
violating phases are not much constrained by current experiment, so we have the freedom 
of taking them quite large. 
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There are more aspects of WIMPy baryogenesis that we have not touched here, but 
we beheve are worth investigating, such as indirect detection signals, and the possibihty 
that the Hght stable state n, into which the exotic quark tp decays, could constitute the 
extra degree of radiation at BBN time. There is also more work related to model-building 
that should be done. For example one could think about UV completions of the models 
presented here, and it would be interesting to explore other discrete symmetry groups that 
would work with the mechanism. 
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A The complete list of dimension six operators 

Given the particle content in Table 1, we want to write down all the possible dimension 
six operators. Our Lagrangian has the form 



In this paper we set the scale A to 10 TeV. The reader should keep in mind that varying 
the couplings Aj's is equivalent to varying the scale A, the only measurable quantity in the 
EFT being the combination Aj/A. 

The couplings Aj's are complex in principle. As we show in Appendix D, we definitely 
need some of them to be complex in order to have the CP violation necessary for the 
generation of a baryon asymmetry. 

The strategy for writing down all possible operators is: 

• write down all the operators consistent with the symmetries, 

• use Fierz identities to reduce to a linearly independent basis. 
Operators linear in Tp and u 



These operators are relevant for the annihilation of DM, proceeding through an s-channel 
(first line) or through a t-channel (second line). 




(A.l) 



+XliX^^){Xu) + Xl{X^^)iXu) 



(A.2) 



+h.c. 
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Operators bilinear in u 

A?(Xn)(Xtu'f) + Xl{Xu)iX^u^) + h.c. (A.3) 

These contribute to the DM annihilation into a pair of quarks and to the tree-level direct 
detection cross section. 

Operators bilinear in tp and u 

A^(V'V)(nn) + Xloi'ilJu)iij^u^) + Xl^i'^^^){uu) + A?2(V^'^u'^)(V^u) + h.c. (A.4) 

The operators 9 and 11 contribute to the washout cross section. 
Operators bilinear in tp 

A?3(XX)(V^V5) + A?4(Xtxt)(V'V^) 
+A?5(X^)(X^) + Xl,{X^^){Xi;) + Xl,{X^^){X^) + A?8(Xt^t)(j^^) + A?9(Xtv^t)(xV;) 

+h.c. 
(A.5) 

These contribute to the DM annihilation only when < . If one wants to successfully 
generate a baryon asymmetry through the annihilation of DM, these operators need to be 
suppressed with the respect to the ones linear in ip and u. 

Operator for the decay of ip 

Xlo{'ipd){dn) + h.c. (A.6) 

Depending on the charge assignment of n, see Table 1, this operator either conserves baryon 
number or violates it. 

B The role of the discrete symmetry 

Discrete symmetries, such as i?-parity in SUSY or KK-parity in Universal Extra Dimen- 
sions, just to mention two popular examples, are a generic feature of models with DM 
candidates. In this appendix we show that, considering the Abelian discrete group Z„, we 
need at least n = 4 for our models. 

Let us start with a generic Z„ and let us assign the following charges 

Qx = exp ( —qx j , = exp ( —q^j j , = exp ( —qn j ■ (B.l) 

To avoid proton decay we charge all the SM quarks, but keep all the leptons neutral. 
This way a baryon cannot decay into mesons and/or leptons. Then we have the following 
requirements: 

• dark matter, X, has to be stable. This implies that qx 7^ 0; 
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• we want to avoid decays of the exotic quark, tp, into SM quarks only, in order not to 
spoil the asymmetry. Given that we have charged the SM quarks, we can avoid such 
dangerous decays by keeping ip neutral, = 0; 

• to generate the asymmetry, we need DM to annihilate both into ip + u and their 
conjugates. This requires two operators, {XX){il)u) and {XX){ip^u'^). 

The last requirement imposes the following conditions on the charges 

2gx + fe = O(modn), (B.2) 
2gx - Qu = (mod n), (B.3) 

from which we find either qx = or qx = n/A. As we said, for DM stability, X has to be 
charged, so the first solution is not acceptable. We conclude that 

qx = n/4, qu = n/2. (B.4) 

Given that the g^'s must be integers, this implies n = ik, with k an integer. This proves 
that the discrete group has to be Z4fc, the smallest one being Z4. 

The solution (B.4) has an important implication: the discrete charge of X is complex. 
This means that a Majorana DM does not work here. Therefore, the simple considerations 
that we have outlined force the dark matter fermion to be Dirac in our models. Note 
that we do not expect this to be a generic feature of WIMPy baryogenesis models: it is 
conceivable that other non-Abelian discrete symmetries can stabilize DM without requiring 
it being Dirac. 

The Z4 symmetry also prevents any coannihilation process between ip and x- 
C Boltzmann equations 

The evolution of the DM, ^ and baryon asymmetry number densities in the expanding 
Universe is governed by a set of Boltzmann equations. Introducing the rescaled inverse 
temperature z = rn^/T and the comoving number densities 1^ = n^{z)/s{z), s{z) being 
the entropy density, we can write: 
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where H{z) is the Hubble expansion rate. The thermahy averaged interaction rates are 
defined for a, i j ■v^ kl scattering by 

= r^s^sK, (^) al]{s), (C.5) 

J Sinf V / 

were the boundary is s^f = max|(mj (m/j + m;)^}, and a is the reduced cross 

sections are given in terms of the Mandelstam variable by 

and the integration limits are given for example in the PDG [63]. 

The Boltzmann equations are quadratic in the DM density Ydm = Yx + Y-^, but we 
only expand them to first order in the asymmetries e and Yaq = Ya — Ya, with a = u,d, ip, 
as these are expected to be small. Up to this approximation and given the list of operators 
above, these equations are the most general ones one can write. They reflect the fact that 
while DM annihilates through many channels, only a part of them {XX ipu + h.c.) is 
CP-violating (CPV) and contributes to the generation of SM and exotic baryons asymme- 
tries. The CPV interactions are proportional to the u-quark number violating operators 
C'i...6, and we provide an approximate formula for in eq. (4.1). The CP-conserving 

annihilations involve the operators 07^8 a-nd C'i3...i9- As the SM baryon-number is con- 
served by the latter processes, they do not contribute to the generation of the BAU. In 
that respect, they constitute a competitive effect that reduces the WIMPy baryogenesis 
efficiency, so they need to be suppressed. We henceforth neglect them. Further, the num- 
ber asymmetries undergo washout processes. The mixed washout, (operators C'i...6), 
mix DM with u-ip, whereas the pure washout, (operators Og and On), involves only 
u and V'- We plot in figure 1 the rates for various limiting cases. 

In equations C.3 and C.4, we introduced the decays of ■0 into two down-quarks, 
parametrized by the rate 7^). This decay provides the only source of d-quark number 
violation. Let us note that the decays of ip inject high-energy down-quarks in the thermal 
bath. To avoid spoiling Big-Bang nucleosynthesis predictions (see, e.g. [64]), these decays 
have to be fast enough, with T[ip ddn) < 1 s. This places constraints on A„ that are 
rather weak. In order to simplify the Boltzmann equations we can ask that ■0 be in thermal 
equilibrium up to the freeze-out of washout. Estimating decay rates by 

^ ddn) - , (C.7) 

we then obtain 

r„ . > .,.wo) ^ A„ > s . X (^) - (^) {^) . 

(C.8) 

The effect of these fast decays are drastic for the abundance: every ip produced in X 
annihilation immediately decays into two anti-down quarks, enforcing An^ ~ 0, while the 
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Figure 12. Diagrams for the annihilation process XX — )■ utp. The diagrams in the first 
two Unes involve s-channel-annihilation operators, while in the last line they involve t-channel- 
annihilation operators. The cross in the last diagram represents a mass insertion, . The annihi- 
lation processes XX — >■ i2^^^ are obtained from similar diagrams, with the replacements A3 — > A|, 
A2 XI, Ae 4-> A5 and Ag -> Xg. 



total baryon asymmetry yielded is approximately An^ + An^ ~ 3 An^. At the end, the 
Boltzmann equations simplify to: 



zsiz)Hiz)'-^ = -2jZ^iz) f (1^) -ij , (C.9) 



zs{z)Hiz)'^ = (^(|^) - ij 

V-'DM 



(C.IO) 



We numerically solve this set of equations. 

D Calculation of the CP asymmetry 

The physical CP asymmetry is defined as 

^ j{XX uil)) + -i{XX uip) + j{XX ntV't) + ry(^xx ut^t) ' ^ ' ' 

where the 7's denote the thermally averaged interaction rates, as defined in the previous 
appendix. Eq. (D.l) can be calculated in our model, at leading order, from the interference 
between tree-level and one-loop diagrams, as shown in figure 12. 
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When we sum the squared matrix elements at the numerator and at the denominator 
we find the following: 



2^|A^NUMp = X |2m2(Xi2+X34) + (s-2m2)(Xi3+X24) 

t — m| 1 

+m^m^{l25 - I35 - 2^16 H 2"(^16 + ^46 2:56)) > , (D.2) 

s - J 

I-^DENP = ^{(^ - 2"^x)(^ - + I^2P + IA3P + IA4P) 



+(t - m2)(t - - m2 )(|A5p + lAeP) 

-4m^(s - m^)(7^i4 + 7^23) - 2m^m^{t - m^)(7^26 + '^45) 

-2m^m^{s + t-w?^- m^)(7^l5 + T^se)} , (D-3) 

with Xjj = Im ^A? A| Ag^^ and TZij = Re ^A? A^^^ A few comments are in order at this 
stage. First, the cross sections at the denominator are computed at tree level. Second, 
the combinations of the couplings that appear in the imaginary parts at the numerator are 
invariant under rephasing of the fields X, tp and n, as one would expect for a physically 
meaningful result. Third, the factor of (s — m^)^ at the numerator is a consequence of the 
fact that the particles in the loop go on-shell. 

Next we need to integrate over the phase space to obtain the cross sections, then we 
need to perform the integrals defined in eq. (C.5). The result we find for e is quite lengthy, 
but, remarkably, analytic. It depends on the couplings A1...6, Ag, on the masses and m^, 
as well as on the temperature via z = m^jT . It is instructive to look at some limiting cases 
rather than at the general full expression. Setting some equalities among the couplings, as 
in Section 2, and taking the low temperature limit (z — )■ cx)), we find 



where as before x = m^/{2m^) and 

Awo = |Awo|e*^ = A9. (D.5) 

With the assumed equalities among the couplings we lose the rephasing invariance of the 
result. This is not a big deal, given that we have shown above that the general result is 
rephasing invariant, but it leads us to make a further assumption, which is to take the 
couplings Ap, At real, and keep only A wo complex. This expression is the same as the one 
we would get by defining e as a ratio of cross sections instead of thermal averages. This is 
no surprise, given that we have taken the zero temperature limit. We would like to remind 
the reader that in the literature e is often defined as a ratio of widths or cross sections, 
instead of thermal averages. The former is an approximation which corresponds to the 
zero-temperature limit of the latter, and is often a good approximation [65, 66]. 
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